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Abstract. Multivariate Poisson distributions have numerous applications. 
Fast computation of these distributions, holding constant a fixed set of linear 
combinations of these variables, has been explored by Sontag and Zeilberger 
|SZ10| . This elaborates on their work. 



1. Introduction 

Set N = Z> . Let X\,... ,X n be independent Poisson random variables with 
Xj ~ Poisson(Xi). Let A G Mat mx „(N). Define new random variables Y\, . . . , Y m 
by taking linear combinations of the Xf. 

(!)"(!) 

Then the joint p.d.f for Yi, . . . , Y m is given: 



F Y (&i, ■ • ■ ,b m ) = P(Yi = 6i, . . . ,Y m = b m ) 



n g -(Ai+A 2 H hA„) 



^ /,-,! /■■>! /,■„' 

(fe!,...,fe„)eN" 

(6i,...,6 m ) T =A(fei,...,fc„) T 

Sontag and Zeilberger demonstrated that the multivariate generating function 
ofY: 

Gy(zi,...,z m ) = ^2 Fy{h,---,b m ) z^ 1 ■■■z%? 

(6i,...,6 m )€N m 

may be simplified [SZIO] (cf. |Burf2j V. 

(n m 
3 = 1 »=1 

Certain conditions can be placed on the matrix A to guarantee that the joint 
p.d.f. Fy can be expressed as the sum over a single index (and thus practicably 
computed with Wilf-Zeilberger recurrences |WZ92[ IAZ06 ). In particular, Sontag 
and Zeilberger established that A £ Mat2 X n({0, 1}) is a sufficient criterion for 
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reduction to a single- index sum |SZ10j . This paper will explore another sufficient 
criterion on the matrix A to guarantee that the joint p.d.f. can be expressed as the 
sum over a single index. 

2. Smith Normal Form 

Let A £ Mat mx „(N). Due to the nature of the problem, it may be assumed 
without loss that the rows of A are linearly independent. In particular, this imposes 
m < n. 

Any matrix A G Mat mx „(Z) can be decomposed into its Smith normal form 
|Jac09| pp. 181-185]: i.e.: there exist invertible (over Z) matrices P e Mat mxm (Z) 
and Q £ Mat nx „(Z) such that: 



V 



\ 



/ 



PAQ 



where r is the rank of A, d\, . . . , d r ^=0, and d\ \ ■ ■ ■ \ d r . These di, 1 < i < r are 
unique up to sign-change and are called elementary divisors. For 1 < i < r, let 
denote the g.c.d. of the i-rowed minors of A. Then d\ — ±Ai, c?2 = iAjAj" 1 , . . . , 
d r = ±A r .A"i 1 . 

Theorem. If the matrix A £ Mat mx „(N) has rank r — m — n— 1 and elementary 
divisors di = ±1 for all 1 < i < r, then Fy^ki, . . . , k m ) can be written as the sum 
over a single index. 



Proof. Note that: 

Fy(bi, . . .,b m ) = 



A^ 1 X k 2 2 



A 



-(Ai+A 2 H hA„) 



fci! fc 2 ! fc«! 



E 

(k 1 ,...,k n )eN rl 
(b 1 ,...,b m ) T =A(k 1 .,...,k n ) T 

By hypothesis, there exist invertible matrices P € Mat„_i x „_i(Z) and Q G 
Mat nx „(Z) such that: 

(l„-i|0» x i) =PAQ 

where l n -i denotes the {n — 1) X [n — 1) identity matrix. Then (pi, ... , b m ) T = 
A(fci, . . . , k n ) T may be rewritten: 



(l)i-i |0„ x i) Q 



'fci' 



\k n i 



so that the first n — 1 entries of the vector Q _1 (fci, . . . , k m ) T agree with the (n — 
1) x 1 vector P(6i, . . . , b m ) T . In particular, by introducing j £ N as the index of 
summation it follows that the full set of solutions to A(fci, . . . , k n ) T = (b\, . . . , b m ) T 
for (fci, . . . , fc„) g N n is given: 



Q 



V o 



p 

1 / 



( h \ 

b n -i 
V 3 J 
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Substituting (ki, . . . , k n ) as given above yields an expression for Fy(bi, . . . , b m ) as 
the summation over a single index j. □ 

Example. Let: 

A _ ( 1 1 1 
JK I 2 1 1 

Then: 

So that A(fci, . . . , k n ) T — (b\, . . . , b m ) T has solutions: 



6)-oa)-(8?i)(5)-( 



6i-6 2 +2j 
J 

b 2 -2j 



So that: 



\fel \fc 2 \fe3 

Oi,fc 2 ,fc 3 )eN 3 

(6 1 ,6 2 ) T =A(fc 1 ,fc 2 ,/ C3 )' r 



A 2 A 3 p -(Ai+A 2 +A 3 ) 

4- (6i 

Example. Let: 



S ( 6l _" &2 + 2 j)!j!(6 2 -2i)! e 



Then: 



A / 1 3 2 2 \ 
A = 5 16 12 17 
V 3 16 21 56 / 

/l \ „.„ /I 00W1 3 2 2 \ 

10 0= PAQ = 1 5 16 12 17 
VOOIO/ \0 1/ \3 16 21 56 / 

So that A(fci, . . . , fc 4 ) T = . . . , & 3 ) T has solutions: 



144 -31 4 15 
-69 15 -2 -5 
32 -7 1 -1 
1 



'fei\ /6i\ / 144 -31 4 15 \ /bi\ / 1446i -31& 2 +4b 3 + 15j 

fe 2 1 - O I b 2 | _ [ -69 15 -2-5 U _ -696i+156 2 -26 3 -5j 

fc 3 — I 6s _ I 32 -7 1 -1 J I 63 _ I 326!-76 2 +6 3 -j 

fc 4 / \ 3 / V 1 / \ J / V j 



So that: 



\fel \fe 2 \fe3 \fc 4 

(fei,fe 2 ,fc3,fe 4 )eN 4 

(&l,6 2 ,6 3 ) T =A(fci,fc 2 ,fc 3 ,fc 4 ) T 
^1446 1 -316 2 +46 3 +15j^-696i+156 2 - 26 3 - 5j^326i-76 2 +6 3 -j^ e -(Ai+A2+A 3 +A 4 ) 

~ ^ (1446i - 316 2 + 46 3 + 15j)!(-696i + 156 2 - 26 3 - 5j)!(326i - 7b 2 + b 3 - Jjljl 

3. Special Case: A is invertible 

Proposition 1. If the matrix A e Mat„ xn (Z> ) is invertible, then i*V(&) can be 
expressed explicitly. 



Proof. Note that: 



fceMat„ xl (Z> ) 
Ak=b 



\ki \k 2 \k n 
Fy(b)= T _2_ . . . _™_ e ~( A i 

Yl J ^ fci! fc 2 ! fcj 



i+A 2 +-+A„) 
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But this sum is taken over the singleton: 

{k e Mat„ x i(Z>o) : Ak = b} = {k £ Mat„ x i : k = A^b) 

□ 



Example. Let: 



. / 1 5 3 

A = 2 10 5 

18 

T 



Then det(A) = 1 and for any b = (b%, 62, ^3) & Mat3 X i(Z>o), the unique solution 
to Ak = b for k G Mat nxl (Z> ) is given: 

kt\ / 75 -37 -5 \ / bi \ /75bi-376 2 -5fc 3 

fc 2 = A b = -16 8 1 62 = -16f)i+86 2 +6 3 

k 3 J \ 2 -1 / \b 3 J \ 2&i-b 2 

Thus the explicit probability density function for Y is given: 



^756i-376 2 -5fc 3 ^-166i+8fc 2 +6 3 ^2fci-h 2e -(A 1 +A 2 +A 3 ) 

(75b! - 376 2 - 56 3 )!(-166i + 8b 2 + 6 3 )!(2&i - 62)! 
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